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ABSTRACT 

20 July 201 1 We present the first direct and unbiased measurement of the evolution of the dust mass func- 

tion of galaxies over the past 5 billion years of cosmic history using data from the Science 
Demonstration Phase of the He rsche I- AlTL AS. The sample consists of galaxies selected at 
250/im which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. 
Dust masses are calculated using both a single temperature grey-body model for the spec- 
tral energy distribution and also using a model with multiple temperature components. The 
dust temperature for either model shows no trend with redshift. Splitting the sample into bins 
of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At 
z = 0.4 — 0.5, massive galaxies had dust masses about five times larger than in the local 
Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the 
optical depth derived from fitting the UV-sub-mm data with an energy balance model was 
also higher This increase in the dust content of massive galaxies at high redshift is difficult to 
explain using standard dust evolution models and requires a rapid gas consumption timescale 
together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or 
combinations of all three. This evolution in dust mass is likely to be associated with a change 
in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier 
cosmic epochs. 
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1 INTRODUCTION 

The evolution of the dust content of galaxies is an important and 
poorly understood topic. Dust is responsible for obscuring the UV 
and optical light from galaxies and thus introduces biases into our 
measures of galaxy properties based on their stellar light (Driver et 
al. 2007). The energy absorbed by dust is re-emitted at longer infra- 
red and sub-millimetre (sub-mm) wavelengths, providing a means 
of recovering the stolen starlight. Dust emission is often used as an 
indicator of the current star formation rate in galaxies - although 
this calibration makes the assumption that young, massive stars are 
the main source of heating for the dust and that the majority of 
the UV photons from the young stars are absorbed and re-radiated 
by dust (Kennicutt et al. 1998, 2009; Calzetti et al. 2007). Many 
surveys of dust emission from 24-850/im (Saunders et al. 1990; 
Blain et al. 1999; Le Floc'h et al. 2005; Gruppioni et al. 2010; Dye 
et al. 2010; Bales et al. 2010) have noted the very strong evolution 
present in these bands and this is usually ascribed to a decrease 
in the star formation rate density over the past 8 billion years of 
cosmic history (z ~ 1: Madau et al. 1996, Hopkins 2004). The 
interpretation of this evolution is complicated by the fact that the 
dust luminosity of a galaxy is a function of both the dust content 
and the temperature of the dust. It is pertinent to now ask the ques- 
tion "What drives the evolution in the FIR luminosity density?", is 
it an increase in dust heating (due to enhanced star formation activ- 
ity) or an increase the dust content of galaxies (due to their higher 
gas content in the past) - or both? 

Dust is thought to be produced by both low-intermediate mass 
AGB stars (Gehrz 1989; Ferrarotti & Gail 2006; Sargent et al. 
2010) and by massive stars when they explode as supemovae at 
the end of their short lives (Rho et al. 2008; Dunne et al. 2009; 
Barlow et al. 2010). Thus, the dust mass in a galaxy should be 
related to its current and past star formation history. Dust is also 
destroyed through astration and via supernovae shocks (Jones et al. 
1994), and may also reform through accretion in both the dense and 
diffuse ISM (Zhukovska et al. 2008; Inoue 2003; Tielens 1998). 
The life cycle of dust is thus a complicated process which many 
have attempted to model (Morgan & Edmunds 2003; Dwek et al. 
1998; Calura et al. 2008, Gomez et al. 2010; Gall, Anderson & 
Hjorth 2011) and yet the basic statistic describing the dust content 
of galaxies - the dust mass function (DMF) - is not well determined. 

The first attempts to measure the dust mass function were 
made by Dunne et al. (2000; hereafter DOO) and Duime & Bales 
(2001 ; hereafter DEOl) as part of the SLUGS survey using a sample 
of IRAS bright galaxies observed with SCUBA at 450 and 850/im . 
Vlahakis, Dunne & Bales (2005; hereafter VDB05) improved on 
this by adding an optically selected sample with sub-mm observa- 
tions. These combined studies, however, comprised less than 200 
objects - none of which were selected on the basis of their dust 
mass. These studies were also at very low-z and did not allow for a 
determination of evolution. A high-z dust mass function was esti- 
mated by Dunne, Bales & Bdmunds (2003; hereafter DEE03) using 
data from deep sub-mm surveys. This showed considerable evolu- 
tion with galaxies at the high mass end requiring an order of magni- 
tude more dust at z ^ 2.5 compared to today (for pure luminosity 
evolution), though with generous caveats due to the difficulties in 
making this measurement. Finally, Bales et al. (2009) used BLAST 
data from 250-500/Ltm and also concluded that there was strong 
evolution in the dust mass function between z — — 1 but were 
also limited by small number statistics and confusion in the BLAST 
data due to their large beam size. 

In this paper, we present the first direct measurement of the 
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space density of galaxies as a function of dust mass out to 2: — 0.5. 
Our sample is an order of magnitude larger than previous studies, 
and is the first which is near 'dust mass' selected. We then use 
this sample to study the evolution of dust mass in galaxies over the 
past ~ 5 billion years of cosmic history in conjunction with the 
elementary dust evolution model of Edmunds (2001). 

The new sample which allows us to study the dust mass func- 
tion in this way comes from the //eric7ie/-Astrophysical Terahertz 
Large Area Survey (H-ATLAS; Bales et al., 2010), which is the 
largest open-time key project currently being carried out with the 
Herschel Space Observatory (Pilbratt et al., 2010). H-ATLAS will 
survey in excess of 550 deg"^ in five bands centered on 100, 160, 
250, 350 and 500/im, using the PACS (Poglitsch et al., 2010) and 
SPIRE instruments (Griffin et al., 2010). The observations consist 
of two scans in parallel mode reaching 5a point source sensitivi- 
ties of 132, 126, 32, 36 and 45 mjy in the 100, 160, 250, 350 & 
500/im bands respectively, with beam sizes of approximately 9", 
13", 18", 25" and 35". The SPIRE and PACS map-making are de- 
scribed in the papers by Pascale et al. (201 1) and Ibar et al. (2010), 
while the catalogues are described in Rigby et al. (2011). One of 
the primary aims of the Herschel-AThAS is to obtain the first un- 
biased survey of the local Universe at sub-mm wavelengths, and 
as a result was designed to overlap with existing large optical and 
infrared surveys. These Science Demonstration Phase (SDP) ob- 
servations are centered on the 9'' field of the Galaxy And Mass 
Assembly (GAMA; Driver et al. 201 1) survey. The SDP field cov- 
ers 14.4 sq. deg and comprises approximately one thirtieth of the 
eventual full H-ATLAS sky coverage. 

In section [2] we describe the sample that we have chosen to 
use for this analysis and the completeness corrections required. In 
section [3] we describe how we have derived luminosities and dust 
masses from the Herschel data, while in section |4l we present the 
dust mass function and evaluate its evolution. Section |6] compares 
the DMF to models of dust evolution in order to explain the ori- 
gin of the strong evolution. Throughout we use a cosmology with 

= 0.27, = 0.73 and = 71kms~^Mpc"\ 



2 SAMPLE DEFINITIONS 

The sub-mm catalogue used in this work is based on the > 5a at 
250/xm catalogue from Rigby et al. (2011), which contains 6610 
sources. The 250/im fluxes of sources selected in this way have 
been shown to be unaffected by flux boosting, see Rigby et al. 
(2011) for a thorough description. Sources from this catalogue are 
matched to optical counterparts from SDSS DR7 (Abazajian et al. 
2009) down to a limiting magnitude of r-modelmag =22.4 using 
a Likelihood Ratio (LR) technique (e.g. Sutherland & Saunders 
1992). The method is described in detail in Smith et al. (2011). 
Briefly, each optical galaxy within 10"of a 250/im source is as- 
signed a reliability, 7?, which is the probability that it is truly as- 
sociated with the 250/im emission. This method accounts for the 
possibility that true IDs are below the optical flux limit, the po- 
sitional uncertainties of both samples, and deals with sharing the 
likelihoods when there are multiple counterparts. For our study we 
have used a reliability cut of R ^ 0.8 as this ensures a low con- 
tamination rate (< 5 percent) which leaves 2423 250/im sources 
with reliable counterparts. The LR method tells us that ~ 3800 
counterparts should be present in the SDSS catalogue, however we 
can only unambiguously associate around 64 percent of these. Our 
sample is thus low in contamination but incomplete (we will deal 
specifically with the incompleteness of the ID process in the next 
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Figure 1. Templates for three galaxies showing the range of optical fluxes 
expected for galaxies which are at the SPIRE flux limit of 5250 = 32 mJy 
at z = 0.5; the limit of our study. The templates are for M82 (a typical 
starburst), a fferac/ie/-ATLAS template derived from our survey data by 
Smith et al. in prep and Arp 220, a highly obscured local ULIRG. The 
SDSS limit of r = 22.4 is shown as a horizontal dotted line and even a 
galaxy as obscured as Arp 220 is still visible as an ID to our optical limit 
at z = 0.5. The yellow shape represent the SDSS-r band filter which was 
used to compute the optical flux 

section). A further cut was made to this sample to remove any stars 
or unresolved objects, this was done using a star-galaxy separation 
technique based on optical/IR colour and size, similar to that used 
by Baldry et al. (2010). Only six objects in the final reliable ID cat- 
alogue have 'stellar or QSO IDs' and so required removal. We also 
removed the five sources which were identified as being lensed by 
Negrello et al. (2010). 

We then used the GAMA database (Driver et al. 201 1) to ob- 
tain spectroscopic redshifts for as many of the sources as possible 
(GAMA target selection is based on SDSS so no further match- 
ing is required). These are supplemented by public redshifts from 
SDSS DR7 (Abazajian et al., 2009), 2SLAQ-LRG (Cannon et al., 
2006), 2SLAQ-QS0 (Groom et al., 2009) and 6dFGRS (Jones et 
al., 2009). Where spectroscopic redshifts were not available we 
used photometric redshifts which were produced for H-ATLAS us- 
ing SDSS and UKIDSS-LAS (Lawrence et al. 2007) data and the 
ANNz method (Collister & Lahav 2004). This method is described 
fully in Smith etal. (2011). 

Section 12.11 shows that we can quantify the statistical com- 
pleteness of the IDs out to 2 = 0.5 and we choose this as the 
redshift limit of the current study. The total number of sources in 
the final sample is 1867 with 1095 spectroscopic redshifts. With 
this sample, the number of expected false IDs (summing 1~ R, see 
Smith et al. 201 1) is 60 (or 3.2 percent). 



2.1 Completeness corrections 

There are three sources of incompleteness in this current sample. 

(i) Sub-mm Catalogue Incompleteness (Cs): This is due to the 
250/im flux limit of the survey and the efficacy of the source ex- 
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Figure 2. Left: SDSS r-modelmag as a function of 250/.imflux. There is no strong correlation apart from at the brightest fluxes. Only 4 galaxies lie within 
0.4 mag of the flux limit used for IDs (r < 22.4) at 2 < 0.5 and so we consider that optical incompleteness is not a serious problem for this sample. Right: 
r-mag versus redshift for all sources in GAMA-9 (pale blue squares) and SPIRE IDs with R > 0.8 (black triangles). Herschel sources tend to be larger mass 
optical galaxies and so the SDSS flux-limit does not affect our ability to optically identify H-ATLAS sources until z ~ 0.5. Note that the light panel uses the 
brighter limit of r < 19 appropriate for the GAMA redshift survey. 



Table 1. The percentage completeness of our reliable ID catalogue as a 
function of redshift, as taken from Smith et al. (201 1). The correction factor 
used in the luminosity function is denoted by ■ 



z 


Completeness (%) 




0.0-0.1 


93.2 


1.07 


0.1-0.2 


83.2 


1.20 


0.2-0.3 


74.2 


1.35 


0.3-0.4 


55.6 


1.80 


0.4-0.5 


53.1 


1.88 



traction process. The catalogue number density completeness has 
been estimated through simulations and presented by Rigby et al. 
(201 1). Apart from the very small range of flux near to the limit, at 
32 — 34 mjy the catalogue is > 80 percent complete. Correction 
factors are applied to each source in turn based on its flux following 
Tables 1 and 2 in Rigby et al. (2011). The largest correction is in 
the flux range 32-32.7 mJy and is a factor 2.17, this applies to 124 
sources out of a total of 1867 at z < 0.5. 

(ii) ID Incompleteness (Cz): The LR method measures in an 
empirical way a quantity Qo, which is the fraction of SPIRE 
sources with counterparts above the flux limit in the optical sur- 
vey. However, it is not possible to unambiguously identify all these 
counterparts with > 80 percent confidence due to positional un- 
certainties, close secondaries and the random probability of finding 
a background source within that search radius. Smith et al. (2011) 
have estimated a completeness for reliable IDs as a function of red- 
shift. This allows us to make a statistical number density correction 
in redshift slices for the sources which should have a counterpart 
above the SDSS limit in that redshift slice, but which do not have 
R ^ 0.8. This correction is applied to each source and is listed in 
Table[T] The ID incompleteness is a function of redshift (not unex- 
pectedly) with corrections of a factor ^ 2 needing to be applied in 
the highest redshift bins. 

(iii) Optical catalogue incompleteness (C ): This correction is 
required because the SDSS catalogue from which we made the 



identifications is itself incomplete as we approach the optical flux 
limit of r = 22.4. We ascertained the completeness using the back- 
ground source catalogue used in the ID analysis of Smith et al. 
(2011), containing all sources which passed the star-galaxy sepa- 
ration at r-modelmag < 22.4 in the primary SDSS DR7 catalogue 
in a region of ~ 35 degrees centered on the SDP field. We fit- 
ted a linear slope to the logarithmic number counts in the range 
r = 19 — 21.5 and extrapolated this to fainter magnitudes. We then 
used the difference between observed and expected number counts 
to estimate completeness. The results are presented in Table[2]and 
show that completeness is above 80 percent to r = 21.8, falling to 
50 percent by r = 22.2. By restricting our analysis to z < 0.5 we 
keep 97 percent of the sources below r ~ 22 and so in the range 
of acceptable completeness. It is possible, in principle, for there to 
be some form of optical incompleteness in the sample which is not 
corrected for with the above prescription, e.g. a population of ob- 
jects which begin to appear at high redshifts in the H-ATLAS sam- 
ple but which are not well represented in SDSS. Such a population 
could conceivably consist of very obscured star-bursts. To test our 
susceptibility to this, we estimate the SDSS r magnitude of a highly 
obscured galaxy with an SED like that of Arp 220 (Av = 15) at 
our 250^m flux limit at the redshift limit of 2: — 0.5 and find that it 
would still be detected in our sample. Figure [T] shows three differ- 
ent SED templates normalised to S250 ~ 32 mJy at z = 0.5: M82, 
an H-ATLAS based template appropriate for sources at z — 0.5 
from Smith et al. in prep, and Arp 220. All templates less obscured 
than Arp 220 are easily visible at our optical flux limit. We will 
therefore proceed on the assumption that no such new populations 
exist below the optical limit in our highest z bins. 

Figure|2ji plots r-mag as a function of 250/im flux. A galaxy with 
5*250 below ~ lOOmJy can have a wide range of optical magnitude 
(r-mag = 16.5-22.0), and while optical magnitude is a strong func- 
tion of redshift this is not the case for the sub-mm flux. Figure [2J5 
shows r-mag as a function of redshift for all galaxies in the GAMA 
9hr (Driver et al. 201 1) spectroscopic sample (cyan), as well as the 
reliable SPIRE IDs (black). This shows a lack of Herschel sources 
at the fainter magnitudes at low redshifts (i.e. the lowest absolute 
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Table 2. The percentage completeness as a function of r magnitude for 
the catalogue used to make the identifications to H-ATLAS sources. The 
connection factor used in the luminosity function is denoted by C'r ■ 



r mag Completeness (%) Cr 



21.6 


91.1 


1.10 


21.7 


87.6 


1.14 


21.8 


82.8 


1.21 


21.9 


77.7 


1.29 


22.0 


70.5 


1.42 


22.1 


61.6 


1.62 


22.2 


52.5 


1.90 


22.3 


42.8 


2.33 


22.4 


17.0 


5.88 



magnitudes or stellar masses)Q It appears that H-ATLAS is less 
sensitive to low stellar mass galaxies than the SDSS (due to them 
having lower dust masses) and so only at high-z does the r-band 
limit preclude the identification of Herschel sources. 

3 DUST MASS AND LUMINOSITY 

The Herschel fluxes are translated into monochromatic rest-frame 
250/im luminosities following 

L250 = 4nD^ (1 + z) S250 K (1) 

where I/250 is in WHz~^, D is the co-moving distance, 525o is the 
observed flux density at 250^m and K is the K-correction which is 
given by: 

where Vohs is the observed frequency at 250/im, Vohs(i+7.) is the 
rest-frame frequency and Tiso and (3 are the temperature and emis- 
sivity index describing the global SED shape. 

In order to derive the values required for the K-correction, a 
simple grey body SED of the form S (x B{v, T) was fitted 
to the PACS and SPIRE fluxes as described in Dye et al. (2010), 
with a fixed dust emissivity index of /3 = 1.5 and a temperature 
range of 10-50 K. Where insufficient data points are available for 
the fit (300/1867), the median temperature of 26 K from the galax- 
ies which could be fitted was used. With only SPIRE data on the 
Rayleigh- Jeans side of the SED (as is the case for most sources), 
only the combination of /3 and Tiso is well constrained, with the 
two parameters being inversely correlated by the fit; good fits are 
obtained with /3 = 1.5 — 2.0. These simple grey body fits can be 
performed for the majority of sources and are accurate at repre- 
senting the flux between rest-frame 250-166/im (relevant for our 
redshift range) and so are suitable for applying the K-correction. 

A dust mass can also be calculated from the observed 
250/im flux density and the grey body temperature as: 



1 The limit of GAMA is r ~ 19 which is brighter than the SDSS limit 
used for H-ATLAS IDs (r ~ 22.4). 



M..^ ^-f/^ + -)^ (3) 

K250 t3\V2Zl3, Jiso j 

where K250 is the dust mass absorption coefficient which we 
take to be equal to 0.89 kg~^ at 250/im (equivalent to scaling 
K850 = 0.077 m^kg"\ as used by DOO, James et al. 2002, da 
Cunha, Chariot & Elbaz 2008 with a /3 = 2). It also lies within the 
range of values found for the diffuse ISM in the Milky Way and 
other nearby galaxies (Boulanger et al. 1996; Sodroski et al. 1997; 
Bianchi et al. 1999; Planck Collaboration 2011a). The dust mass 
via Eq.[3]scales as Md oc y~2.4 ^j- ^ ^ Q for temperatures around 
20 K; changing the temperature from 20-30 K results in a reduction 
in mass by a factor 2.6. At z = 0.5 this dependence is steeper since 
the peak of the dust emission is shifted to longer wavelengths so 
the observed frame is even further from the Rayleigh- Jeans regime. 
Changing from /3 = 1.5 to 2.0 reduces the temperatures by ^ 3K 
and increases the dust masses by ~ 30 — 50 percent. 

This isothermal dust mass estimate can be biased low as it is 
now well established that dust exists at a range of temperatures in 
galaxies. Only dust in close proximity to sources of heating (e.g. 
star forming regions) will be warm enough to emit at A ^ lOO/im 
but this small fraction of dust (by mass) can strongly influence the 
temperature of the isothermal fits. The bulk of the ISM (and there- 
fore the dust) resides in the diffuse phase which is heated by the 
interstellar radiation field to a cooler temperature typically in the 
range of 15-20 K (Helou et al. 1986; DEOl and references within; 
Popescu et al. 2002; VDE05; Draine et al. 2007; Willmer et al. 
2009; Bendo et al. 2010; Boselli et al. 2010; Kramer et al. 2010; 
Bernard et al. 2010; Planck Collaboration 2011b). For more accu- 
rate dust mass estimates we require the mass-weighted temperature 
of the dust emitting at 250/im which requires fitting a model with 
multiple (at least two) temperature components. This is not to say 
that the FIR fluxes for most of the H-ATLAS galaxies are not fitted 
adequately by the single temperature model; an isothermal model 
and a more realistic multi-temperature model are often degenerate 
in their ability to describe the SED shape with a limited number 
of data points. To illustrate this, we show in Fig [3] an example of 
isothermal and 2-component SED fits to a H-ATLAS source with a 
well sampled SED. Although the 2-component fit is formally bet- 
ter, there is nothing to choose between them as descriptions of the 
fluxes of the H-ATLAS source between 60-500/im . DEOI studied 
this issue for a sample of SLUGS galaxies with 450/im detections, 
and concluded that the best overall description of that sample was 
a two-temperature model with /3 = 2 and a cold component tem- 
perature of ~ 20K. 

To deal with the cold dust component, we now introduce a 
more sophisticated SED model which includes dust in several phys- 
ically motivated components, following the prescription of Chariot 
& Fall (2000). The results of this fitting are presented and described 
in detail in Smith et al. in prep, and outlined here in brief. This sim- 
ple, but empirically motivated, SED model fits broadband photom- 
etry from the UV-sub-mm to estimate a wide variety of parameters 
(da Cunha et al. 2008 - hereafter DCE08; da Cunha et al. 2010a). 
The method uses libraries of optical and infrared models (25,000 
optical and 50,000 infrared) and fits those optical-IR combinations 
which satisfy an energy balance criteria to the data. The optical 
libraries have stochastic star formation histories and the stellar out- 
puts are computed using the latest version of the Bruzual & Char- 
lot (2003) population synthesis code (Chariot & Bruzual in prep) 
libraries and a Chabrier (2003) Galactic-disc Initial Mass Function 
(IMF). 
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Figure 3. Top row: Isothermal and 2-component SEDs for an H-ATLAS sources with a well sampled SED. Redshifts and fitted pai'ameters are shown in each 
panel. For the isothermal fits T and /3 were free to vary while for the 2-component fits /3 was fixed to be 2. The parameter Nc/N^ is the ratio of cold/warm 
mass. 



The dust mass in this model is computed from the sum of 
the masses in various temperature components contributing to the 
SED, including cool dust in the diffuse ISM, warm dust in birth 
clouds, hot dust (transiently heated small grains emitting in the 
mid-IR) and PAHs. In the fits to H-ATLAS sources (and SINGS 
galaxies; DCE08) around 90 percent of the dust mass is in the cold 
diffuse ISM component and this is also the best constrained com- 
ponent due to the better sampling of the FIR and sub-mm part of 
the SED with Herschel. Many other studies also find that the cold 
dust component dominates the overall mass, and so it is the most 
important one to constrain when measuring the dust mass function 
(e.g. DEOl; VDE05; Draine et al. 2007; Willmer 2009; Liu et al. 
2010). The priors used by DCE08 for the temperatures of the grains 
in equilibrium are 30-60 K for the warm component and 15-25 K 
for the cold component. These values agree well with temperatures 
measured for local galaxies (Braine et al. 1997; Alton et al. 1998; 
Hippelein et al. 2003; Popescu et al. 2002; Meijerink et al. 2005, 
DEOl, V05, Stevens et al. 2005, Stickel et al. 2007; Draine et al. 
2007; Willmer et al. 2009; Planck Collaboration 2011b) and also 
with temperatures measured from stacking optically selected galax- 
ies with the same stellar mass and redshift range as our sample into 
Herschel-KIhAS maps (Bourne et al. in prep). The value of k used 
in the DCE08 model is (by design) comparable with that used in 
the isothermal fits here. 

The prior space of the parameters is sampled by fitting to sev- 
eral million optical-FIR model combinations and returns a proba- 
bility density function (PDF) for the dust mass and other param- 
eters (e.g. dust temperature, stellar mass, dust luminosity, optical 
depth and star formation rate) from which the median and 68 per- 
cent confidence percentiles are taken as the estimate of the quantity 
and its error. 

This model was fitted to the 60 percent of the galaxies in our 
sample for which useful optical and NIR data were available from 
GAMA. We fitted only to galaxies which have matched aperture 
photometry in r-defined apertures as this best represents the total 
flux of the galaxy in each band as described in Hill et al. (2011); 
Driver et al. (2011). The distribution of sources with and without 
these SED fits as a function of redshift is shown in Figure|4] Those 
without fits dominate only in the highest redshift bin, from z — 
0.4-0.5. 

The errors on the dust mass range from ±0.05 — 0.27 dex 
and this error budget includes all uncertainties in the fitting from 
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Figure 4. The distribution of sources with DCE()8 SED fits as a function of 
redshift (red), those without fits are shown in blue (dashed). 



flux errors to changes in temperature and contribution of the vari- 
ous dust components. Some typical SED fits and PDFs for the dust 
mass and cold temperature parameters are shown in Figure m The 
dust mass is generally a well constrained parameter of these model 
fits; the PDF is narrower when more IR wavelengths are available 
and so the cold temperature is then better constrained. 

A comparison of the isothermal dust masses ( Miso ) and the 
full SED based masses ( Mscd ) is shown in Figure |6^ and there 
is generally poor agreement between the two, with the scatter of a 
factor 2-3 related to the difference between the temperature of the 
isothermal fit and that of the DCE08 model fit. This sensitivity is 
because at 250/im we are near to the peak of the black body func- 
tion for the cold temperatures appropriate to the bulk of the dust 
mass (15-20 K). At longer sub-mm wavelengths (such as 850/im ), 
this temperature sensitivity is less severe, but the choice of dust 
temperature used when estimating masses at rest wavelengths close 
to those of Herschel is clearly important. 

For sources which have insufficient UV-sub-mm data to use 
the DCE08 model, we need to extrapolate dust masses by com- 
paring 1/250 with Msed for those sources which do have fits. The 
relationship is linear, with some scatter introduced by the range in 
dust temperature for the cold ISM component (Figure |6j3): 

log Mscd = log L250 - 16.47 (4) 



Evolution of dust mass 7 




Wavelength / microns ]oq{UjM^) if/K 




Wavelength / microns \aq(Ug/Ma) if 





Figure 5. SED fits and probability distribution functions for dust mass and diffuse ISM dust temperature for a range of H- ATLAS sources. The black curve on 
the SED plot is the total attenuated starlight and re-radiated dust emission. Blue curve is the unattenuated starlight. Green is the attenuated starlight and red is 
the dust emission. The red squares show the observed photometry and en'ors or upper hmits. The Hmit to the dust mass accuracy is our ability to determine 
the cold temperature, which is better constrained when there are more FIR data points available. The best constraints on the dust mass are ~ 0.05 dex and the 
worst are ~ 0.27 dex. 



Equation |4] is used to convert L250 to dust mass in cases where 
the full SED could not be fitted (747 sources out of 1867). The re- 
lationship between M^cd and the cold temperature of the diffuse 
ISM (which dominates the dust mass in these galaxies) is similar to 
that in Eqn. [S] since the DCE08 model fits the sum of grey-bodies 
at different temperatures to the photometry. The colder the temper- 
ature fitted, the higher the dust mass will be for a given I/250 ■ This 
is clearly demonstrated in Fig. |6] (bottom). In using this relation- 
ship for sources without SED fits we are making the assumption 
that they will also fall on this relationship and that there is no sys- 
tematic trend in those sources without fits (e.g. if only the highest 
redshift or most/least luminous sources did not have fits). Since the 
galaxies without fits span a range of redshift and luminosity and as 
we also find no correlation of temperature with L250 for our sample 
(see Section |42l and Figure [T4t. there should be no bias introduced 
in using the relationship in Figure [6J3 to estimate masses for those 
sources without fits. 

The scatter in this figure is influenced by our choice of prior 
for the temperature of the diffuse dust component (15-25 K). Al- 



lowing a wider prior will broaden the scatter if a significant number 
of sources are best fitted by hotter or colder temperatures. This is- 
sue is explored further by Smith et al. in prep who conclude that for 
a sample of galaxies with well constrained cold temperatures this 
prior encompasses ^80 percent of the population. Further study of 
the temperatures of the populations requires a larger sample with 
good 5 band FIR/sub-mm photometry, which will be possible with 
the next Phase of H- ATLAS data comprising 10 times the area of 
SDP. The difficulty in using a wider prior is that when we lack full 
coverage of the FIR/sub-mm SED (as is the case where we have 
only limits at PACS wavelengths and 500/im ) there is only a weak 
constraint on the cold temperature (aTc ~ 2.5 — 3 K). This can 
place a galaxy with a real temperature of 15 K down at 12 K, and 
produces quite a bias in the fitted dust temperature (since at 12 K 
the mass is very sensitive to temperature). The model will fit ar- 
bitrarily high masses of very cold dust since this contributes very 
little to the overall energy balance. Our choice to restrict the tem- 
perature prior to the parameter space which is preferred by obser- 
vations of nearby galaxies, and by those galaxies well sampled in 
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Figure 6. Top: Comparison of dust masses from the DCE08 model ( A/ged ) 
versus the dust mass obtained by using the isothermal grey body fit ( Miso ) 
using Eq. [3] Points are colour-coded by the isothermal temperature. The 
one-to-one line is shown in black. There is a large difference between the 
two mass estimates which is a strong function of fitted isothermal temper- 
ature. Bottom: Comparison of Mgcd to L250 showing a reasonably tight 
and linear correlation. The best fit relationship (Eq. |4) is over-plotted. The 
scatter in this relationship is driven by the diffuse ISM dust temperature, 
which is used to colour-code the points. 



H- ATLAS, potentially means we underestimate the masses of some 
cooler sources, but we would prefer to be conservative at this point. 



4 THE DUST MASS FUNCTION 
4.1 Estimators 

To calculate the dust mass function we use the method of Page & 
Carrera (2000; hereafter PCOO) who describe a method to estimate 
binned luminosity functions that is less biased than the 1/Vmax 
method (Schmidt 1968). To begin with, we produce measurements 
of the 250/im luminosity function since this is more directly re- 
lated to the flux measurements from Herschel and enables us to 
discuss the method without the added complication of translating 
250/xm luminosity into dust mass. The PCOO estimator is given by 



EAf ^ ^ 



az 



(5) 



where Cs, Cz Cr are the completeness corrections for each 
object as described in Section[2jT]and the sum is over all galaxies in 
a given slice of redshift and luminosity bin. Lmax and Lmin are the 
maximum and minimum luminosities of the bin. Zmin is the mini- 



mum redshift of the slice and 2max(L) is the maximum redshift to 
which an object with luminosity, L, can be observed given the flux 
limit and K-correction, or the redshift slice maximum, whichever 
is the smaller. The PCOO method has the advantage of properly cal- 
culating the available volume for each L — z bin and, in particu- 
lar, it does not overestimate the volume for objects near to the flux 
limit. This prevents the artificial turn-down produced by 1 / Vmax 
in the first luminosity bin of each redshift slice. We compare to the 
1/Vmax estimator in Figure |7] and confirm that the 1/Vmax esti- 
mate of the 250/xm luminosity function suffers from the bias noted 
by PCOO due to slicing in redshift bins. 

In the PCOO formalism described above, the accessible volume 
is not calculated individually for each source (as for 1/Vmax) but 
is instead calculated for each bin in the L — 2 plane using a global 
K-correction. However, we know that each object in our sample 
has a different K-correction because they have different grey body 
SED fits. We therefore modified the estimator such that the acces- 
sible volume for a given L ~ z bin is calculated for each galaxy in 
that bin in turn using its grey body SED fit to generate its limiting 
redshift Zmax,i ~ z{Li, Smin, Td) across the bin. These individual 
contributions are then summed within the bin such that: 



^ ^ rijmax rz 
'=1 JL„i„ Jz, 



Cs Cz Cr 



Sj^dz dL 

az 



(6) 



Note that this is not the same as reverting to the 1 /Vmax esti- 
mator as we are still calculating the volume available for each L~z 
bin, however we are now being more precise about the shape of the 
limiting curve for each source based on its individual SED. This is 
clear from the difference in the LF calculated this way, as shown in 
Figure |7lb) compared to the PCOO and 1/Vmax methods shown in 
Figure |7{ a). This change affected the highest redshift bins most as 
expected. 

In this case, the error on the space density is given by 



(70 



(7) 



where is the individual (j) contribution of a galaxy to a par- 
ticular redshift and luminosity bin, and the sum is over all galaxies 
in that bin. The error bars in Figure|7]show these errors. 

This 250/im luminosity function differs slightly from that pre- 
sented in Dye et al. (2010) in that the ID sample has since been up- 
dated to include extra redshifts (1867 compared to 1688) and also 
to remove stars, for which there were 130 contaminating the pre- 
vious sampl^ While Dye et al. (2010) did attempt to correct for 
incompleteness in the optical IDs of the sub-mm sample, we are 
now able to extend this to correct for incompleteness as a function 
of redshift, r-mag and sub-mm flux which was not previously pos- 
sible. The results are, however, comparable in that strong evolution 
in the 250/im LF is evident out to z ^ 0.4. There is then seemingly 
a halt, with little evolution between z = 0.4 and z = 0.5. This is 
still consistent with Dye et al. (2010) within the error bars of both 
estimators. 

We suspect that this behaviour in the highest redshift bins is a 
result of a bias in the ANNz photo-z we are using. Figure[8^ shows 
a comparison between spectroscopic and photometric redshifts of 



Due to using an earlier version of the LR estimate which combined stars 



and galaxies together 
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Figure 7. Top: 250/im luminosity functions calculated via the 1/Vmax method (open triangles /dashed) and PCOO method (solid circles and lines) in five 
redshift sUces of Az = 0.1 out to 2 = 0.5. Colours denote the redshifts as: black (0 < z < 0.1), red (0.1 < 2 < 0.2), green (0.2 < 2 < 0.3), blue 
(0.3 < 2 < 0.4) and cyan (0.4 < 2 < 0.5). The bias in 1/Vmax in the lowest luminosity bin in each redshift slice is apparent from the turn-down in this 
bin. Bottom: 250^m luminosity function calculated using the modified PCOO estimator which includes an individual K-correction for each object in an L — 2 
bin. Using individual K-corrections has a more significant effect in the highest redshift slices as expected. 
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Figure 8. Top: Comparison of spectroscopic versus photometric redshifts 
for galaxies in tlie H- ATLAS SDR There is a bias in the photo-z at redshifts 
greater than ~ 0.3 where photo-z tend to underestimate the true redshift. 
A 1-to-l correlation is shown by the solid line. Bottom: The fraction of 
photo-z used in the luminosity function as a function of redshift. Photo-z 
start to dominate the LF at the same redshift where the photo-z bias begins. 

H-ATLAS sources in the SDP region. There is a bias above z ~ 
0.3 — 0.35 where the photometric redshifts tend to underestimate 
the true redshift (see Fleuren et al. in prep). This issue is further 
exacerbated by the fact that this is also the redshift at which the LF 
becomes dominated by photo-z (Fig.HJ)). 

There is another potential bias in the highest-z shce due to 
the optical flux limit approaching the main body of galaxies in the 
sample. While we correct for the incompleteness in space density 
due to the r-band limit, we are not able to deal with any accompa- 
nying bias which might allow only those galaxies with lower dust- 
to-stellar mass ratios into the sample at the highest redshifts (see 
Fig EH and Section [6] for more discussion). Greater depth in opti- 
cal/IR ancillary data will be required to test the continuing evolu- 
tion of the luminosity and dust mass functions beyond 2 = 0.5 and 
this will soon be available with VISTA- VIKING and other deep op- 
tical imaging for the H-ATLAS regions from VST-KIDS, INT and 
CFHT. 

Having demonstrated that our modified version of the PCOO 
estimator produces sensible results on the 250/im luminosity func- 
tion, we now turn to the estimate of the dust mass function (DMF). 
We again use Eqn |6] however we now sum all galaxies in a bin 



of the Md — z plane. We use the ratio of Md to L250 to estimate 
the Lniax and Lmin for each galaxy, which is required to compute 
the individual K-correction. The results for both single temperature 
masses and SED based masses are shown in Figure|9] 

Both estimates of the dust mass function show a similar evolu- 
tionary trend as the 250^m LF, with the same apparent slow down 
at higher redshift which we believe may be related to issues with 
the photo-z. The evolution is present whichever estimate of the dust 
mass is used, however, we will continue the discussion using the 
DMF from the DCE08 SED based dust masses (Fig.|9)i) as we be- 
lieve that this is the best possible estimate at this time. 

The dust mass function also shows a down-turn in some red- 
shift slices at the low mass end. We do not believe that this rep- 
resents a true dearth of low mass sources at higher redshifts but 
rather reflects the more complex selection function in dust mass 
compared to I/250 . While there is a strong linear relationship be- 
tween our dust mass and L250 (Figure [6J5) there is still scatter on 
this relationship due to the variation in the temperature of the cold 
dust in the ISM. At fixed L250 warm galaxies will have smaller dust 
masses than cooler ones, which leads to a sort of 'Eddington' bias 
in the dust masses. At the limiting 1/250 for a given redshift bin we 
are not as complete as we think for low dust masses, since we can 
only detect galaxies with low dust masses if their dust is warmer 
than average. This in turn leads to the apparent drop in space den- 
sity. In the two highest redshift bins, the fraction of sources without 
SED fits increases and so the dust masses are then directly propor- 
tional to the 250/im luminosity. To improve on this, we would need 
to use a bi-variate dust mass/ L250 approach for which the current 
data are insufficient, however this analysis will be possible with the 
complete H-ATLAS data-set. 

4.2 Dust temperatures and evolution of the LF 

An important ingredient in our estimate of the dust mass is the dust 
temperature. In order to interpret the increase in sub-mm luminos- 
ity as an increase in the dust content of galaxies, we have to be 
wary of potential biases in our measurements of the dust temper- 
ature. The dust temperature is most accurately constrained when 
PACS data, which span the peak of the dust SED, are available in 
conjunction with the longer sub-mm data from SPIRE, which also 
constrains the cold temperature. For the SDP field, the PACS data is 
shallow and this results in PACS detections for only 262 galaxies. 
The fraction of PACS detections as a function of redshift are 41, 21, 
8, 7, 4 percent respectively from the lowest to highest redshift bin. 
A comparison of the temperatures from fits with PACS detections 
and without in the lowest redshift bin (where PACS samples a rep- 
resentative fraction of the population) is shown in Figure [TO] Fits 
without PACS detections are included only if the 350/im flux was 
greater than 3ct in addition to the > 5ct 250^m flux. The left panel 
shows the cold dust temperature from the DCE08 fits and it can 
be seen that PACS sources have a range of cold ISM temperatures, 
however, those without PACS detections tend to have mostly cooler 
dust in their ISM. Smith et al. in prep show that the DCE08 fit- 
ting does tend to underestimate the cold temperature slightly when 
PACS data are removed from a fit, but this effect is of order 1-2 K 
and does not fully account for the difference in these distributions. 
When we consider instead the peak temperature of the SED, that 
given by the isothermal fit (right panel in Fig.[TO), we see a differ- 
ent trend. Now the PACS detections are found only at the higher 
end of the temperature range while those without PACS detections 
span a wider range of temperature. There is no bias when the PACS 
data are removed from the fits (the temperatures vary randomly by 
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Figure 9. Dust mass functions using the modified PCOO estimator calculated in 5 redshift slices of Az = 0.1. Top SED based dust masses. Bottom Isothermal 
dust masses. The relation between dust mass and L250 has some scatter due to variation in the temperature of the cold ISM dust, which results in down-turns 
in the lowest mass bins in each redshift slice. The broader error on A/j acts to convolve the true DMF with a Gaussian of width approximately 0.2 dex. 
Schechter functions are plotted in the top panel with the faint-end slope fixed to that which fits best in the z < 0.1 slice. Parameters for the fits are given in 
Table[3] 
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~ ±3 K when the PACS data are removed). The sources with the 
coldest isothermal temperatures (Tiso < 20K) are not detected by 
PACS, even in the lowest redshift bin, as they either do not contain 
enough warm dust, or are not massive enough, to be detected in our 
shallow PACS data. We also note that where PACS does not pro- 
vide a > 5(T detection, we use the upper limit in the fitting which 
provides useful constraints on dust temperature for many more H- 
ATLAS sources. 

The trend for PACS to detect the warmer sub-population of 
H- ATLAS sources becomes more pronounced at higher redshift, as 
the galaxies must be intrinsically more and more luminous to be 
detected by PACS. Figure [TT] shows the fitted temperature versus 
redshift for sources which are detected by PACS (black points) and 
those not-detected by PACS but which have at least 2 good quality 
sub-mm points for the fit. If we used only the PACS detected sam- 
ple, we would infer an evolution in dust temperature in this redshift 
range - but this is a selection bias due to the sensitivity of the PACS 
bands to warm dust and the shallow survey limit for PACS. 

We can also look at the dust colour temperatures of the H- 
ATLAS sources in comparison to other samples without the com- 
plications of fitting models. In Figure[T2]we compare the FIR/sub- 
mm colours of the 35 H- ATLAS sources which have 60, 100 and 
500/im detections at ^ 3cr with the colours of SLUGS galax- 
ies from DEOl and VDE05 to see how these sub-mm selected 
sources compare to those selected at 60/imfrom the IRAS BGS 
(Soifer et al. 1989) or in the optical. H- ATLAS fluxes at 60/^mare 
from the IIFSCz catalogue of Wang & Rowan-Robinson (2009), 
1 00/im fluxes are from PACS and 500^mfrom SPIRE (Rigby et 
al. 2011). To allow a comparison between 500/im fluxes from H- 
ATLAS and 450^m fluxes from SLUGS, we reduce the SLUGS 
450/im values by 37 percent using a standard template suitable for 
SLUGS sources from DEOl (approximately oc v"^ at these wave- 
lengths). All of these sources are local. Figure [72l shows that the 
H-ATLAS sources are significantly colder in their colours than the 
warmest end of the IRAS sample; they overlap rather better with 
the optically selected SLUGS sample. This is not surprising given 
our selection at 250//m is more sensitive to the bulk dust mass of 
a galaxy while that at 60/im from IRAS is more sensitive to warm 
dust (either large, warm grains in star forming regions or small tran- 
siently heated grains). We note that, since only a very small number 
(35) of H-ATLAS sources are detected by IRAS, these few sources 
shown in Fig.[T2]are also likely to have 'warmer' colours than the 
overall H-ATLAS sample. This agrees with the findings that PACS 
is sensitive to only the warmer H-ATLAS sources at higher red- 
shifts and the SED fitting results which show that the H-ATLAS 
sources contain relatively cooler dust. 

If the evolution in the 250/im LF were due simply to an in- 
crease in the 'activity' of galaxies of the same dust mass, then we 
should see a corresponding increase in dust temperature with red- 
shift and no evolution in the DMF. To explain the amount of evo- 
lution in the 250pim LF without any increase in dust mass would 
require an increase in the average dust temperature of order a fac- 
tor 2 over the period < z < 0.5. We investigated the relationship 
between both the cold ISM dust temperature from the DCE08 fits 
and the isothermal grey-body temperature with redshift and found 
no trend for either (Figure [Tsj at 2 < 0.5, similar to the results 
from Amblard et al. (2010) and inconsistent with the temperature 
evolution required to explain the increase in the 250/im luminosity 
density. 

The temperature of nearby {z < 0.1) dusty galaxies has 
been shown to be correlated with their IR luminosity (the so-called 
Lir — Td relation; e.g. DOO, Dale et al. 2001). A natural explana- 



tion for this observation might be that a galaxy which has hotter 
dust (for a given mass) will have a larger IR luminosity than a sim- 
ilar mass galaxy with cooler dust. Recent work which extends to 
more sensitive surveys and samples selected at longer wavelengths 
suggests that this does not hold at higher redshifts and that galax- 
ies are in general cooler at a given IR luminosity than previously 
believed (Coppin et al. 2008; Amblard et al. 2010; Rex et al. 2010, 
Symeonidis et al. 2009, Seymour et al. 2010, Smith et al. in prep). 
Symeonidis et al. (201 1) suggest that this is due to more rapid evo- 
lution of "cold" galaxies over the period 0.1 < z < 1 than "warm" 
ones. Recent studies at other wavelengths (70-160/im from Spitzer 
and PACS) seem to support this interpretation, finding that cold 
galaxies are responsible for most of the increase in the IR luminos- 
ity density over the range < z < 0.4 (Seymour et al. 2010; Grup- 
pioni et al. 2010). This is in agreement with the evolution seen in 
H-ATLAS galaxies which are largely comprised of this 'cold' pop- 
ulation. Despite our average luminosity increasing with redshift, 
we see no increase in the average temperature (either isothermal 
or cold ISM temperature) and indeed we also see no correlation of 
either temperature with luminosity (either dust luminosity from the 
DCE08 model or L250 ) for this sample (see FigurefT?!. 

To summarise, while we are subject to uncertainties in our 
ability to derive the dust masses and the exact scale of any evo- 
lution, we are nevertheless confident that: 

• The evolution in the 250/xmLF out to 2 = 0.5 cannot be 
driven by dust temperature increases; there must be some evolu- 
tion in the mass of dust as well. 

• The H-ATLAS sources at 2 ^ 0.5 are colder than previous 
samples based on IRAS data and therefore most of the evolution 
at low redshift is driven by an increase in the luminosity or space 
density of such cooler galaxies. 

• H-ATLAS sources show no trend of increase in dust tempera- 
ture with either redshift or luminosity at 2 < 0.5 

4.3 Comparison to low redshift dust mass functions 

We can compare the lowest redshift bin in the DMF (0 < 2 < 0.1) 
to previous estimates from the SCUBA Local Universe and Galaxy 
Survey (DOO; VDE05), which used SCUBA to observe samples of 
galaxies selected either at 60^m from the IRAS Bright Galaxy Sam- 
ple (Soifer et al. 1989) or in the B-band from the CfA redshift sur- 
vey (Huchra et al. 1983). The IRAS SLUGS galaxies were mostly 
luminous star-bursts, and in principle this should have produced an 
unbiased estimate of the local dust mass function as long as there 
was no class of galaxy unrepresented in the original IRAS BGS 
sample. However, it was argued in DOO and VDE05 that this selec- 
tion at bright 60^m fluxes quite likely missed cold but dusty galax- 
ies, given the small sample size of ^ 100, thus may have produced 
a DMF which was biased low. The optically selected SLUGS sam- 
ple overcame the dust temperature bias and did indeed show that 
there were very dusty objects which were not represented as a class 
in the IRAS BGS (similarly confirmed by the ISO Serendipity Sur- 
vey: Stickel et al. 2007). The directly measured DMF presented 
by VDE05 suffered from small number statistics, and instead V05 
followed the work of Serjeant & Harrison (2005) in extrapolating 
the IRAS PSCz (Saunders et al. 2000) out to longer wavelengths 
(850^m ) using the empirical colour-colour relations derived from 
the combination of IRAS and optically selected SLUGS galaxies. 
This set of 850//m estimates for all IRAS PSCz sources was then 
converted to a dust mass assuming a temperature of 20 K (the aver- 
age cold component temperature found by DEOl and VDE05) and 
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Figure 10. Left: Cold ISM temperature from the DCE08 fits for the lowest redshift bin which has 41 percent PACS detections. The sources with PACS 
detections are shown by the red dashed line while those without PACS detections but which do have a 350^m flux above 3cr in addition to the 5cf 250/im point 
are shown in black. Right: Same but for the isothermal temperature 




Figure 11. Left: Cold ISM temperature from the DCE08 fits versus redshift for sources with PACS detections (black filled) and which have 350/xm fluxes 
above 3(t in additions to the 5a 250/xmflux (red open). Right: Same but for the isothermal temperature. Here there is a correlation between Tiso and redshift 
for the PACS detections (r = 0.4). 
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Figure 12. Colour plots for the 35 H-ATLAS galaxies with detections at 60, 100 and 500/im compared to those for SLUGS sources detected at 450/imfTOm 
an IRAS 60/xm selected sample (IRAS) and an optically selected sample (OS). The SLUGS points have had their 450^m fluxes adjusted downward by 37 
percent to make them equivalent to 500/^m 
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Figure 13. Left: The temperature of the cold interstellar dust component as a function of redshift z. Only sources with either a PACS detection or a 35()/^m flux 
above 3a, in addition to the 250fj,m flux, are plotted. Mean values and l-cr errors on the mean are shown as black points. The data points in magenta show the 
full distribution of the temperatures. The large en'or bai' in the top right shows the average 68 percent confidence range on the temperature for an individual fit. 
Right: The isothermal temperature estimated from a grey body fit versus redshift, same coding as before. The line plotted shows the evolution in temperature 
required in order to explain the evolution in the 250/im LF without any increase in the dust masses. Neither method for estimating the dust temperature shows 
any evolution with redshift. 
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Figure 14. Left: Cold dust temperature and L250 showing no correlation. The points are colour coded by redshift. Right: Same as (left) but for the isothermal 
dust temperature 



a mass opacity coefficient of ksso = 0.077 kg ^ . From this set 
of masses they then produced an estimate of the DMF. 

The DMFs are compared in Figure [TS] where the black solid 
line and points are from H- ATLAS at z < 0.1, the blue dot-dash 
line and filled triangles is the SLUGS IRAS directly-measured DMF 
(DOO) and the red dashed line and open triangles is the DMF based 
on the extrapolation of the IRAS PSCz by VDE05. In this figure, the 
H- ATLAS DMF has been corrected for the known under-density of 
the GAMA-9hr field relative to SDSS as required when comparing 
to an all-sky measurement such as SLUGS or IRAS PSCz. This 
correction is a factor of 1.4 (Driver et al. 2011). The SLUGS DMFs 
have been corrected to the cosmology used in this paper, however 
these corrections are small at low-z. 

It is remarkable that despite the considerable differences in 
sample size, area and selection wavelength, the SLUGS estimate 
from VDE05 based on extrapolating the IRAS PSCz gives a very 
good agreement to our measure. This implies that there is not a 
significant population of objects in the PSCz sample, or the H- 



ATLAS sample which is not represented by the combined opti- 
cal and 60/im selected SLUGS samples (which comprised only 200 
objects). Note that had VDE05 used the IRAS data alone to measure 
dust masses, the results would be extremely different. It is only that 
SLUGS allowed an empirical statistical translation between IRAS 
colours and sub-mm flux and from there, assumed a mass-weighted 
cold temperature for the bulk of the dust that they were able to ob- 
tain such a good measure of the DMF. 

The original direct measure of the DMF from the bright IRAS 
SLUGS sample (blue line in Fig [Tsj DOO) dramatically under- 
estimates the dust content in the local Universe (this was also noted 
by VDE05 once the optically selected sample was included). The 
dust masses were derived for those objects in an identical way to 
the VDE05 DMF (and very similar to our current method which has 
an average measured cold temperature of between 15-19 K), how- 
ever the IRAS BGS simply missed objects which were dusty but did 
not have enough warm dust to make it above the 60/im selection. 
Herschel is able to select sources based on their total dust content, 
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Figure 15. Comparison of the local dust mass functions al z < 0.1 from 
H- ATLAS (black solid line and points) along with estimates from SLUGS. 
Blue dot-dash line and solid triangles - directly measured DMF from IRAS 
SLUGS sample (DOO, DEOl). Red dashed line and open triangles - extrapo- 
lated DMF from IRAS PSCz using sub-mm colours from the optical SLUGS 
sample (VDE05). The H-ATLAS points have been corrected for the factor 
1 .4 under-density in the GAMA-9hr field for this redshift range compared 
to SDSS at large. 



rather than simply the small fraction of dust heated to > 30 K. 
Herschel samples are therefore likely to contain a far wider range 
of galaxies in various states of activity, so long as they have enough 
material in their ISM. 



4.4 Evolution of the dust mass function 

For illustration, we now fit Schechter functions (Schechter 1976) 
to the dust mass functions in each redshift slice. Only in the first 
redshift bin do we fit to the faint end slope a, for other redshift 
bins we keep this parameter fixed at the value which best fits the 
lowest redshift bin (a = —1.01) to avoid the incompleteness prob- 
lem mentioned above with the lowest mass bins at high redshift. 
The best-fitting parameters for the slope a, characteristic mass Af^ 
and normalisation are given in Table[3] where the errors are cal- 
culated from the 68 percent confidence interval from the con- 
tours. For the lowest redshift bin, we include errors which reflect 
the marginalisation over the un-plotted parameter The contours 
for and are shown in Figure [T6l 

There is a strong evolution in the characteristic dust mass Af J 
with redshift, from Af* = 3.8 x lO'^ Mq at 2 < 0.1 to = 
3.0 X 10** Mq at z = 0.4 - 0.5. There is seemingly a decline in (f>* 
over the same redshift range, from 0.0059 — 0.0018 Mpc"'^ dex^^ 
(however this could also be due to sample incompleteness which 
is not corrected for despite our best attempts). The drop in <j}* and 
increase in A/J are correlated (see Fig|16l), and therefore we caution 
against using the increase in the fitted alone as a measure of the 
dust mass evolution. If we keep (f)* fixed at 0.005 Mpc~^ dex~^ 
(which is the average of that for the first two redshift bins) then the 

of the highest redshift bins decreases to 1.8 x 10* M0 giving 
an evolution in A/J over the range z = — 0.5 of a factor ~ 5 
rather than ^ 8 as is the case if the normalisation is allowed to 
drop. 

We calculate the dust mass density in redshift slices using 
Eqn.m 
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Figure 16. confidence intervals at 68, 90, 99 percent for and 
with fixed a for the five redshift bins. This shows the clear evolution of 
over the interval < z < 0.5. 



Pd = r(2 + a) Ma (p* 



(8) 



This assumes that we can extrapolate the Schechter function be- 
yond the range over which it has been directly measured. Given the 
low value of a used (~ — 1) the resulting integral is convergent and 
so whether we extrapolate or not has negligible effect on the result- 
ing mass density values. The values for pd are listed in Table|3]and 
shown as a function of redshift in Figure[T7] There is clearly evolu- 
tion in the cosmic dust mass density out to z ~ 0.4 of a factor ~ 3 
which can be described by the relationship pd oc (1 + z)'^'"'. In the 
highest redshift bin the dust mass density appears to drop (despite 
the increase in M^], but we again caution that this may be due to 
incompleteness/photo-z bias in the final redshift bin. This measure 
of the dust mass density at low redshift can be compared to that 
made by Driver et al. (2007). They used the optical B-band disk 
luminosity density from the Millennium Galaxy Catalogue scaled 
by a fixed dust mass-to-light ratio from Tuffs et al. (2004). Their 
quoted value for the dust density is pd = 3.8 ± 1.2 x 10^ Mpc~^ 
at 2 < 0.1 but this is for a k value from Draine & Li (2001) 
which is lower than that used here by 70 percent. Scaling their re- 
sult to our K, and correcting the density of our lowest redshift bin 
by the factor 1.4 from Driver et al. (2011) (to allow for the under- 
density of the GAMA-9hr field relative to SDSS at z < 0.1) we 
have values of pd = 2.2 ± 0.7 x 10^ Mpc~^ (optical based) and 
pd = 1.4 ± 0.2 X 10^ Mpc"^ (DMF) which are in rather good 
agreement given the very different ways in which these estimates 
have been made. 

We can also calculate the dust mass density parameter Sldust 

from 

- P'' 

Pcrit 



fid 



is the critical density for 

6 



where Pcrit = 1.399 x 10^' Mq Mpc 

h = 0.71. This gives values of fidust ~ 0.7 — 2 x 10"" depend' 
ing on redshift. Fukugita & Peebles (2004) estimated a theoretical 
value of f2dust = 2.5 x 10~® today based on the estimated density 
of cold gas, the metallicity weighted luminosity function of galax- 
ies and a dust to metals ratio of 0.2. This is a little higher than our 
(density corrected) lowest redshift estimate of 1.0 ± 0.14 x 10~* 
but not worryingly so. Menard et al. (2010) also estimate a dust 
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Table 3. The Schechter parameters fitted to the dust mass function 
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The first line of the table is the fit to all three parameters for the lowest redshift bin with associated errors from the 68 percent confidence interval derived 
from the contours. The following entries are where a is fixed to the best-fitting value in the lowest redshift bin. The final entry is the fit to the z ~ 2.5 
DMF from DEE03 corrected to this cosmology and K250. Cos. Var. is the cosmic variance estimated using the calculator from Driver & Robotham (2010). 
Nbin is the number of sources in that redshift bin and Zphot / ztot is the fraction of photometric redshifts in that bin. 
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Figure 17. Integrated dust mass density as a function of redshift for H- 
ATLAS calculated using Eqn[8] The best fitting relationship excluding the 
higher redshift point is over-plotted, which isp£joc(l + 2)^'^. 
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Figure 18. Comparison of the H- ATLAS dust mass function in five redshift 
slices (as in Fig|9) and the high redshift. z ~ 2.5, DMF from D03 (magenta 
dashed line). 



density in the halos of galaxies through a statistical measurement of 
reddening in background quasars when cross-correlated with SDSS 
galaxies. They estimate a dust density of = 2.1 x 10~® for 

a mean redshift of 2 ~ 0.35 and suggest that this is dominated 
by 0.5 L, galaxy halos. Comparing this to our measure of the dust 
within galaxies at the same redshift (^f^l^ = 2 x lO"'^) we see 
that at this redshift there is about the same amount of dust outside 
galaxies in their halos as there is within. We note here that dust in 
the halos of galaxies will be so cold and diffuse that we will not 
be able to detect it in emission with H-ATLAS and so it is not in- 
cluded in our DMF. The decrease in pd at recent times could be 
due to dust being depleted in star formation, destroyed in galaxies 
by shocks or also lost from galaxies (and from our detection) to the 
halos. We will return to this interesting observation in Section[6] 

We can compare the DMF from H-ATLAS to that at even 
higher redshifts, as traced by the 850/im selected SMG popula- 
tion. An estimate of the DMF for these sources at a median redshift 
of z ~ 2.5 was presented in DEE03, using the method. 
In Fig [TS] we show this higher- 2; DMF alongside the H-ATLAS 



data, where the z ~ 2.5 DMF is the magenta solid line with filled 
triangles. The DEE03 higher-z DMF has been scaled to the same 
cosmology and value of «:25o as used here. At z; ~ 2.5, observed 
850/.tm corresponds to rest-frame ^ 250/im and so our lower-z H- 
ATLAS sample and the one at z ^ 2.5 are selected in a broadly 
similar rest-frame band. DEE03 used a dust temperature of 25 K to 
estimate the dust mass, which allowed for some evolution over the 
low-z SLUGS value of 20 K. The z ~ 2.5 sources from DEE03 are 
all ULIRGS and these higher luminosity sources do show enhanced 
dust temperatures in the local Universe (Clements, Dunne & Eales 
2010; da Cunha et al. 2010b). It is also consistent with the cold, 
extended dust and gas component (T = 25 — 30 K) of the highly 
lensed SMG al z — 2.3 (Swinbank et al. 2010; Danielson et al. 
2011) and other lensed sources discovered by Herschel (Negrello 
et al. 2010). If we were to recompute the z ~ 2.5 DMF using a 
temperature of 20 K instead, this would shift the points along the 
dust mass axis by a factor ~ 1.7. 

For either temperature assumption, the z ^ 2.5 DMF is 
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broadly consistent with tlie H- ATLAS DMF in the two highest red- 
shift bins {z = 0.3 — 0.5). The fits to the high-z DMF are shown 
in Table [3] and the dust density at z ~ 2.5 is also consistent with 
that in the z = 0.3 — 0.5 range from H-ATLAS. If true, this im- 
plies that the rapid evolution in dust mass may be confined to the 
most recent 4—6 billion years of cosmic history. Notwithstanding 
the earlier statement that this trend needs to be confirmed with a 
larger sample, dust masses are unlikely to continue rising at this 
pace because the dust masses at very high redshifts (Michalowski 
et al. 2010; Pipino et al. 2010) are not very different to those we see 
here. 

This implies that the evolution in the 250/^m LP is due at least 
in part to a larger interstellar dust content in galaxies in the past 
as compared to today, at least out to z ~ 0.4 (corresponding to a 
look-back time of 4 Gyr). However, an increase in star-formation 
rate is also an important factor as if the dust mass increased at a 
constant SPR we would expect to see a decline in dust temperature 
with redshift. Our observations thus point to an increase in both 
dust mass and star formation activity. If the evolution in the DMP 
is interpreted as pure luminosity (or mass) evolution (as opposed 
to number density evolution), then this corresponds to a factor 4-5 
increase in dust mass at the high mass end over the past 4 Gyr. Since 
dust is strongly correlated to the rest of the mass in the interstellar 
medium (ISM) (particularly the molecular component), this also 
implies a similar increase in the gas masses over this period. In 
contrast, we know that the stellar masses of galaxies do not increase 
with look-back time, showing very little evolution in the mass range 
we are dealing with (predominantly L, or higher) (Pozzetti et al. 
2007; Wang & Jing 2010). The evolution of the DMP is therefore 
telling us something quite profound about the evolution of the dust 
content of galaxies, and by inference, the gas fractions of galaxies 
over this period. 



5 THE DUST CONTENT OF H-ATLAS GALAXIES 

There are two ways in which we can quantify the dust content: 
the amount of light absorbed by dust (or opacity), and the dust- 
to-stellar mass ratio. Both of these are derived from the DCE08 
SED model fits for galaxies which were bright enough (r ^5 20.5) 
that aperture matched photometry was extracted by GAMA (Hill 
et al. 2011). Due to this being shallower than the depth to which 
we can ID the H-ATLAS sources we have to take care not to intro- 
duce selection biases when making these comparisons. Pigure[T9l 
shows r-mag as a function of redshift for the H-ATLAS sources and 
again highlights that H-ATLAS does not detect low stellar mass (or 
low absolute Mr) sources. The panels in Pig.[T9]have colour coded 
points for sources where SED fits were made, and the colours rep- 
resent either the V-band optical depth (top) or the dust-to-stellar 
mass ratio (bottom). At z ~ 0.35 the optical sample which has 
SED fits becomes incomplete, with only the brighter fraction of the 
galaxies having SED fits at a given redshift. This can lead to a low- 
ering of the average optical depth, or dust-to-stellar mass ratio in 
bins aiz > 0.35, since the brighter galaxies (higher stellar masses) 
tend to have lower values of optical depth or dust-to-stellar mass. 
Thus in the following discussion we limit our model comparisons 
to the data with 2 < 0.35. We hope to extend the SED fitting to the 
fainter sources in future work. 

Pirst we plot the amount of optical light obscured by dust: the 
V-band opacity. This is derived from the DCE08 SED model fits, 
and is calculated both in the birth clouds where stars are bom (tv 
from DCE08) and also in the diffuse ISM (fiTy from DCE08). Pig- 




Figure 20. Upper red points: Mean V-band optical depth in the birth clouds 
(from the DCE08 SED fits of Smith et al. in prep) as a function of redshift 
with the best linear fit. Lower black points: V-band optical depth in the ISM 
(fiTv from DCE08). 



ure [20] shows the evolution of both forms of V-band optical depth 
from the model fits, indicating that galaxies are becoming more ob- 
scured back to z ~ 0.4. Choi et al. (2006), Villar et al. (2008) and 
Garn et al. (2010) also find a higher dust attenuation in high redshift 
star forming galaxies. This is sometimes attributed to an increase 
of SPR with look-back time (Gam et al. 2010) and an attendant in- 
crease in dust content rather than to a change in dust properties. It 
is also possible that the apparent increase of optical depth with in- 
creasing redshift is related to the correlation between IR luminosity 
and dust attenuation (Choi et al. 2006), whereby more IR luminous 
galaxies tend to be more obscured. The average IR luminosity of 
our sample increases strongly with redshift (due both to the flux 
limit of the survey and the strong evolution of the LP) and it is cur- 
rently not possible for us to disentangle the effects of redshift from 
those of luminosity since we do not have a large enough sample to 
make cuts in redshift at fixed luminosity. Regardless of which is the 
driver, the observational statement remains that a sub-mm selected 
sample will contain more highly attenuated galaxies at higher red- 
shifts. This is in contrast to some UV selected samples which show 
either no trend with redshift or a decline of attenuation at higher-z, 
due to their selection effects (Burgarella et al. 2007; Xu et al., 2007; 
Buat et al. 2009). This just highlights the obvious - that PIR and 
UV selected samples are composed of quite different objects. 
Our relationships with redshift are as follows: 

birth clouds : rv ~ 3.43z + 1.56 

diffuse ISM : /irv = 1.50z + 0.36 

which implies that the attenuation from the birth clouds is rising 
faster with increasing redshift than that in the diffuse ISM. At 
higher redshifts we are therefore finding that the birth clouds are 
producing a larger fraction of the attenuation in the galaxy than at 
low redshift. We find this trend interesting but further work is re- 
quired to explain and confirm it, firstly ensuring in a larger sample 
that it is not again related to the luminosity (more luminous sources 
also have higher relative attenuation from the birth clouds). Includ- 
ing Balmer line measurements in the DCE08 fits will also better 
constrain the optical depth in the birth clouds. 

Secondly, we can look at dust and stellar mass together us- 



18 L. Dunne et al. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 



Figure 19. r-mag versus redshift for the H-ATLAS sources. Black open circles represent H-ATLAS sources which are too faint for an SED fit using the 
DCE08 model at the cun'ent time, or which were not in the region covered by GAMA photometric catalogues. Coloured points denote the values of either 
V-band optical depth (top) or dust-to-stellar mass ratio (bottom) from the DCE()8 fits. The limit of reasonable completeness in the optical for the SED fits is 
z ~ 0.35. Beyond this redshift, averaged values of optical depth or dust-to-stellar mass ratio will be biased low because only the brightest optical galaxies in 
that redshift bin will have SED fits (and these tend to have less obscuration). 



ing the stellar masses from the DCE08 SED fits. Figure 1271 shows 
the variation of dust and stellar mass with redshift, where the dust 
mass has been scaled up by a factor 178 in order to roughly make 
Md and M, equivalent at the lower boundary at low-z. Magenta 
points show stellar mass, open black squares are the scaled dust 
mass. The stellar mass remains fairly constant with redshift, while 
there is a distinct lack of high dust mass objects in the local Uni- 
verse (as is shown also by the DMF). The dust-to-stellar mass ratio 
as a function of redshift is shown in Figure [23] and discussed in 
more detail in the next section. 



6 MODELLING THE EVOLUTION OF DUST 

In this Section we will attempt to explain the evolution we see in 
the dust content of H-ATLAS sources and in the DMF. We do this 
using a chemical and dust evolution model which traces the yield 
of heavy elements and dust in a galaxy as its gas is converted into 
stars. A full treatment of the evolution of galaxies will be consid- 
ered in Gomez et al. in prep. Here we will consider the elemen- 
tary model of Edmunds (2001; see also Edmunds & Eales 1998) in 
which one assumes that the recycling of gas and dust in the inter- 
stellar medium is instantaneous. Details of the model are given in 
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Figure 21. Stellar mass (magenta) and dust mass scaled by 178 (black open 
squares) versus redshift. The dust mass is scaled to make the dust and stellar 
lower limits approximately coincide at low-z. This illustrates the different 
trends of dust and stellar mass with redshift, with the dust mass evolving 
more rapidly than the stellar mass (as is also evident from the DMF). At 
lower redshifts there are many galaxies with higher stellai' masses than the 
scaled dust mass, while at high redshifts both stellar and dust masses are 
comparable with the same scaling. 

Appendix A, but in brief, a galaxy is considered to be a closed 
box with no loss or addition of gas during its evolution. The evo- 
lution of the galaxy is measured in terms of /, its gas fraction, 
which represents the fraction of the baryonic mass in the form of 
gas. Gas is converted into stars using a star formation prescription 
ipit) — where g is the gas mass and fc is the star formation 

efficiency (inversely proportional to the star formation time-scale). 
We define an effective yield p = p' /a ~ 0.01 where a ~ 0.7 
is the mass fraction of the ISM locked up in stars (EgllOt and p 
is the yield returned from stars for a given initial mass function 
(IMF). We can interpret p as being the true mass fraction of heavy 
elements returned per stellar generation, since some fraction of the 
generated heavy elements is locked up in low mass stars and rem- 
nants. In the first instance, we use the Scalo form of the IMF (Scalo 
1986) for Milky Way evolution (e.g. Calura et al. 2008). The metal 
mass fraction of a galaxy is tracked through p and therefore follows 
metals incorporated into long lived stars and remnants or cycled 
through the ISM where they are available to be made into dust. The 
parameters which determine how many of the available metals are 
in the form of dust relate to the sources of dust in a galaxy and we 
consider three of these: 

(i) Massive stars and SNe: xi is the efficiency of dust conden- 
sation from new heavy elements made in massive star winds or 
supernovae. 

(ii) Low-intermediate mass stars (LIMS): X2 is the efficiency 
of dust condensation from the heavy elements made in the stellar 
winds of stars during their RG/AGB phases. 

(iii) Mantle growth in the ISM: We can also assume that grains 
accrete at a rate proportional to the available metals and dust cores 
in dense interstellar clouds (Edmunds 2001). e is the fraction of 
the ISM dense enough for mantle growth, r/c is the efficiency of 
interstellar depletion in the dense cloud (i.e. if all the metals in the 
dense clouds are accreted onto dust grains then rjc = 1)- 

Morgan & Edmunds (2003) used observations of dust in low- 
intermediate mass stars to show that X2 ~ 0.16 yet theoret- 
ical models following grain growth in stellar atmospheres (e.g. 
Zhukovska et al. 2008) suggest higher values of X2 ~ 0.5. We 



adopt the higher value here, but note that there is some considerable 
uncertainty on X2- For core-collapse supernovae (using theoretical 
models of dust formation e.g. Todini & Ferrara et al. 2001) Morgan 
& Edmunds suggest that xi ~ 0.2; this agrees with the highest 
range of dust masses published for Galactic supernova remnants 
(Dunne et al. 2003; 2009, Morgan et al. 2003; Gomez et al. 2009). 
If core-collapse SNe are not significant producers of dust (e.g. Bar- 
low et al. 2010) or if most of their dust is then destroyed in the 
remnant (Bianchi & Schneider 2007) then this fraction decreases 
to Xi ^ 0-1' making it difficult to explain the dust masses we see 
in our Galaxy or in high-redshift submillimetre bright galaxies with 
stellar sources of dust (e.g. Morgan &. Edmunds 2003; Dwek et al. 
2007; Michalowski et al. 2010). 

For mantles we arbitrarily set e = 0.3 and from interstellar 
depletion levels in our Galaxy and following Edmunds (2001), we 
set r]c ~ 0.7 (that is, we assume that if the clouds are dense, then it 
is likely that the dust grains accrete mantles). In this scenario, the 
dust is formed during the later stages of stellar evolution and uses 
up the available metals in dense clouds. The addition of accretion 
of metals onto grain cores with the parameters described here will 
double the peak dust mass reached by a galaxy. Assuming no de- 
struction of grains, a closed box model and mantle growth gives the 
highest dust mass attainable for galaxies. 

Dust destruction can be added to this elementary model by as- 
suming some fraction 5 of interstellar grains are removed from the 
ISM as a mass ds is forming stars. We use two destruction sce- 
narios: one with a constant destruction rate 5 — 0.3 (Edmunds 
2001) and the second where S is proportional to the Type-II SNe 
rate (which gives a similar result to Dwek's approximation for MW 
IMF; Dwek et al. 2011). We also allow a mantle growth propor- 
tional to SFR since one would expect that the efficiency will de- 
pend on the molecular fraction of the ISM (which in turn is related 
to the SFR; Papadopoulos & Pelupessy 2010). 

Finally, we relax the closed-box assumption and include out- 
flows in the model (Appendix A) since galactic-scale outflows are 
thought to be ubiquitous in galaxies (Menard et al. 2010 made a re- 
markable detection of dust reddening in the halos of galaxies which 
implies at least as much dust is residing in the halos as in the disks). 
Here we test outflows in which enriched gas is lost at a rate propor- 
tional to one and four times the SFR (more powerful outflows are 
unlikely, since in the latter case, the galaxy would only retain ap- 
proximately 20 percent of its initial gas mass). 

6.1 Evolution of Dust to Stellar Mass 

The dust-to-stellar mass ratio of the models discussed here is shown 
in Figure [22] over the life-time of the galaxy as measured by the 
gas fraction, /. The shaded region shows the range of values of 
M^/Mt estimated for the H- ATLAS galaxies, which have a peak 
value of 7 X 10^'^ at z — 0.31 and then decreases as the galaxy 
evolves in time (to lower gas fractions) to 2 x 10""^. This global 
trend is reproduced by the closed box model where dust is con- 
tributed by both massive stars and LIMS, or via mantle growth, 
however the models struggle to produce values of Al^/M, as high 
as observed. We also plot in Figure (22] the variation of M^/M^ 
if low-intermediate mass star-dust is the only stellar contributor to 
the dust budget (xi = 0, X2 = 0.5). It is clear that the LIMS 
dust source cannot reproduce the values of dust/stellar mass seen in 
the H-ATLAS sources alone. Either significantly more dust is con- 
tributed to the ISM via massive stars/SNe than currently inferred, 
or a significant contribution from accretion of mantles in the ISM is 
required (indeed we would need significantly more dust accretion 
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in the ISM than dust produced by LIMS). The simple model also 
suggests that the H-ATLAS galaxies must be gas rich (/ > 0.4) 
in order to have dust-to-stellar mass ratios this high. (Typical gas 
fractions for spiral galaxies today are / ~ 0.1 — 0.2.) 

We can also consider the evolution of dust-to-stellar mass as 
a function of time (Eg. I21t. This is shown in Fig|23ti using dust 
production and yield parameters appropriate for spiral galaxies like 
the Milky Way {p = 0.01, a = 0.7, Xi = 0.1, X2 = 0.5, ery^ = 
0.24, k = 0.25 Gyr~^). We compare the model for two formation 
times of 2 = 0.6 and 2 = 1, where formation time in this model 
can simply mean the time of the last major star formation event. In 
this scenario, we would expect any previous star formation to have 
already pre-enriched the ISM with some metallicity Zi, therefore 
increasing the available metals for grain growth in the ISM. 

From Fig|23b, we see that the MW model does not match the 
variation of dust/stellar mass from H-ATLAS observations even if 
we increase the mantle growth or the amount of dust formed by 
stars, since the increase in dust-to-stellar content with gas fraction 
(as we look back to larger redshifts and earlier times in the evolu- 
tion of the galaxy) is simply not rapid enough. Fig 123b shows the 
same two formation times but now we have tuned the parameters 
to match the data for a formation at 2 = 0.6. In order to do this we 
have to increase the SF efficiency parameter (k = 1.5 Gyr~^) to 
produce a steeper relationship as observed. An increase in k com- 
pared to the MW model is hardly surprising, since these higher 
values are typical of star-forming spirals with initial SFRso of 
ij) ~ 50 M© yr~^ which is in agreement with the observations of 
H-ATLAS sources at higher redshifts. However, increasing k then 
dramatically reduces the actual dust content at any epoch due to re- 
moval of the ISM through the increase in star formation efficiency. 
To explain the high A/d /A/* values for the H-ATLAS sample, we 
would then need to increase the dust condensation efficiencies (i.e. 
the amount of metals which end up in dust) to a minimum of 60 per- 
cent and the effective yield p of heavy elements from stars would 
need increase by at least a factor of two. This is much higher than 
observed condensation efficiencies for LIMS or massive stars/SNe 
although the difference could come from mantle growth. An in- 
crease in the effective yield can only be achieved through the IMF. 
The stellar masses of H-ATLAS galaxies are based on the Chabrier 
IMF (Chabrier 2003), which has 0^0.6 (compared to a ~ 0.7 for 
Scalo). However, to significantly increase the yield from the stellar 
populations, we would require a top-heavy IMF (e.g. Harayama, 
Eisenhauer & Martins 2008). In comparison to the MW-Scalo IMF, 
the effective yield p can increase by a factor of 4 and more material 
is returned to the ISM {a < 0.5). A model with these 'top-heavy' 
parameters is shown in Figure l23b (solid blue), and reproduces the 
H-ATLAS observations without the need for extremely efficient 
mantle growth or higher dust contribution from SNe. A top-heavy 
IMF also frees up more gas and metals in the ISM throughout the 
evolution of the galaxy with time, i.e. / ~ 0.5 at 2 = 0.4 compared 
to the / ~ 0.3 for a Scalo IMF, providing a consistent picture with 
the observed high dust-to-stellar mass ratios and the expected high 
gas fraction for H-ATLAS sources. 

If we assume an earlier formation time, or time since last star 
formation phase, the model cannot reproduce the H-ATLAS ob- 
servations and would require even more extreme values for the 
dust condensation efficiency and/or yield. This suggests a time for 
the last major star formation episode for H-ATLAS galaxies to be 
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Figure 22. Variation of dust-to-stellar mass ratio as a function of gas 
fraction. The shaded box region is the range of values observed for the 
H-ATLAS galaxies. The models are (i) a closed box with no gas enter- 
ing/leaving the system with dust from both supemovae xi = 0.1 and 
LIMS stars X2 = 0.5 (thick solid; black); (ii) with dust from LIMS only 
XI = 0, X2 = 0.5 (thin solid; black); (iii) model (i) now including man- 
tle growth (dot-dashed; black); (iv) A model with mantle growth, where 
the mantle rate is proportional to the SFR (solid; red); (v) and (vi) a model 
which has outflow with gas lost at a rate proportional to one or four times 
the SFR (A/q) (dashed; blue). 



somewhere in the past 5-6 Gyr (which is consistent with the de- 
tailed SED modelling of Rowlands et al. in prep). 

In summary, from this simple model, it is difficult to explain 
the high dust-to-stellar mass ratios in the H-ATLAS data even by 
assuming we are observing these galaxies at their peak dust mass 
unless (i) the fraction of metals incorporated into dust is higher 
(although we would require x > 70 per cent of all metals to be in- 
corporated into dust) or X > 50 per cent with pre-enrichment; (ii) 
The yield is significantly increased via a top heavy IMF. An IMF of 



the form ( 



m 1 oc m 



would increase the yield and hence dust 



mass by a factor of four, easily accounting for the highest Md /M* 
ratios. Such IMFs have been postulated to explain observations of 
high-2 sub-mm galaxies, highly star-forming galaxies in the local 
Universe and galaxies with high molecular gas densities (Baugh 
et al. 2005; Papadolpoulos 2010; Gunawardhana et al. 2011). (iii) 
H-ATLAS galaxies are rapidly consuming their gas following a rel- 
atively recent major episode of star formation (at 2 ^ 0.6). 

6.2 Evolution of the DMF 

We now turn to the evolution of the dust mass itself as evidenced 
from the DMF (Fig|9j which shows an increase in the dust mass of 
the most massive sources of a factor 4—5 in a relatively small time- 
scale (0 < 2 < 0.5, At < 5 Gyr). To show the maximum change 
in dust mass in galaxies in the model, we plot the ratio (R) of dust 
mass at time t to that at the present day, assuming a gas fraction of 
/ ~ 0.1 today (Figure [241. For a closed box model, there is little 
evidence for the dust mass in a given galaxy changing by more than 
a factor of 1.5 in the past compared to its present day value. 

It is clear that including outflows produces a better fit to the 
variation of dust mass observed in the DMF, with the maximum 
change in dust mass approaching the observed change in DMF with 
i? ~ 4 for the extreme outflow model. However, in this case, the 
peak Md /M* is at least an order of magnitude below the observed 
values predicting only 2 x 10~* (see Fig|22t. In this scenario, we 
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Figure 23. Left: The dust-to-stellar mass ratio as a function of redshift. Stellar and dust masses are derived from the SED fits using the models of DCE08 and 
ai'e discussed in detail in Section|3]and Smith et al. (201 lb). Black points show those sources with spectroscopic redshifts, while red points include photometric 
redshifts. Each sample is limited in redshift to the point where the optical flux limit is not biasing the selection to low dust-to-mass ratios. The model lines for 
the dust model (Section. 16. It corresponding to the Milky Way including mantle growth and destruction are over-plotted with formation redshifts of 2 = 0.6 
(dot-dashed) and z = 1 (dotted). A model including pre-enrichment of Zi ~ O.IZq with formation timescale at 2 = 0.6 is also shown (solid; black). 
Right: Same as left including pre-enrichment, but models are now tuned to match the data for the z = 0.6 formation time. With pre-emichment, we require 
Xi = 0.1, X2 = 0.5, p = 0.02, e = 0.9 and SF efficiency fc = 1.5 Gyr~^ to 'fit' the data points (black dot-dashed) or Xi = X2 = £ = 0.5, p = 0.02 (not 
shown). Also shown is a model with mantle growth varying with SFR and a top-heavy IMF described by a = 0.5, p = 0.03 (solid; blue). Adding outflow or 
destruction rates which vary with SFR would make the decline in M^/Mt, more pronounced at lower redshifts (later evolutionary times). 



Closed E=5=0 

Closed Xi=0,X2=0.5 

Closed £=0.3, 6=0.3 

Closed e a SFR 




Gas Fraction, f 

Figure 24. Ratio (R) of dust mass at gas fraction / to that at / = 0.1 
(today). The models are (i) a closed box with no gas entering/leaving the 
system with dust from both supernovae xi =0.1 and LIMS stars X2 = 0.5 
(thick solid; black); (ii) with dust from LIMS only xi = 0, X2 = 0.5 
(thin solid; black); (iii) including mantle growth (dot-dashed; black); (iv) A 
model with mantle growth proportional to the SFR (solid; red), (v) and (vi) 
A model which has outflow with gas lost at a rate proportional to one or 
four times the SFR (A/a) (dotted; blue). It is worth noting that for higher 
returned fraction from stars to the ISM (i.e. a = 0.5), the ratio decreases 
for all models (R < 3 for the extreme outflow). 

would require x > 0.8, erj > 0.8 and p > 0.03. Such high dust 
condensation efficiencies from stellar sources are not observed in 
the MW, and a yield as high as p = 0.03 would again, imply a 
top heavy IMF. For an outflow model with A/a = 1.0, the pa- 
rameters X > 0.6, £77 > 0.3 and p > 0.02 would be required to 
produce the H-ATLAS dust-to-stellar mass ratios, these are more 
reasonable values yet this outflow rate is not sufficient to account 
for the increase in dust mass seen in the DMF (reaching a maxi- 



mum R ~ 1.5; Fig|24t. We believe that outflows must be present 
at some level (Alton, Davies & Bianchi 1999) and the observation 
made earlier that there is as much dust in galaxy halos as there is 
in galaxies themselves is strong circumstantial evidence for some 
outflow activity. Given that there are other ways (e.g. radiation pres- 
sure on grains; Davies et al. 1998) to remove dust from disks, we 
can attempt to derive a rough upper limit for the outflow required 
to produce as much dust in halos at z --^ 0.35 as found by Menard 
et al. (2010). We integrate the dust mass lost from outflows during 
the evolution of the galaxy and compare this to the dust mass in 
the galaxy at z = 0.3 — 0.4 for various values of outflow and star 
formation efficiency k. The results are shown in Table |4] This as- 
sumes no dust destruction in either the halo or the disk, and as such 
is a very simple model. Equality in dust mass inside and outside 
galaxies can be achieved by 2; = 0.3 by having moderate outflow 
< 4 X SFR and 0.25 < k < 1.5Gyr~^. This is not to say that 
all galaxies need have similar evolution; it is quite likely that IT- 
ATLAS sources are more active and dusty and as such may contain 
more dust in their halos than the average SDSS galaxy probed by 
Menard. This simple exercise merely gives some idea of what sort 
of 'average' chemical evolution history is required to reproduce the 
observation. 

We now have a conundrum in that the observed evolution in 
dust mass requires significant outflow of material, however such 
outflow leads to the lowest values of dust-to-stellar mass ratio and 
cannot be reconciled to the observations without extreme alter- 
ations to the condensation efficiencies for dust or the stellar yields. 
Including dust destruction and mantle growth models which vary 
with the SFR alleviates this somewhat since both decrease the dust 
mass more significantly at later times. The change in dust mass over 
the same period compared to the elementary model with constant 
e and S is then more pronounced, but not enough to explain the 
evolution in the DMF. 

One solution to this is if the galaxies with the highest dust 
masses at 2 0.4 — 0.5 are not the progenitors of the H-ATLAS 
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Table 4. t is the age since formation of tlie galaxy at z = 0.6. Outflow 
= 1 and 4 is outflow proportional to 1 and 4 times the star formation rate. 
'Halo/Disk' is the ratio of the integrated dust mass lost in outflow from 
iform to t divided by the dust mass in the galaxy at t. 



z 


t 

(Gyr) 


Outflow = 
fc = 0.25Gyr-l k 
Halo/Disk 


1 

= 1.5Gyr-i 
Halo/Disk 


Outflow = 4 
k = 0.25Gyr-l 
Halo/Disk 


0.5 


0.5 


0.09 


0.42 


0.33 


0.4 


1.0 


0.2 


0.96 


0.73 


0.3 


2.2 


0.4 


3.03 


1.95 


0.2 


3.2 


0.5 


5.47 


3.24 


0.1 


3.5 


0.6 


12.2 


5.13 



sources at « ~ 0.1. We speculate on a scenario where the low 
redshift spiral galaxies {z < 0.15) which do fit the MW model 
in Fig|23b comprise one population and the higher redshift (more 
dusty) objects are a rapidly evolving star-burst population with 
much higher star formation efficiencies (higher k), higher dust con- 
densation efficiencies and/or top-heavy IMFs. The fate of the high 
redshift dusty population is that they rapidly consume their gas (and 
dust) in star formation and by low redshift they are no longer de- 
tected in H-ATLAS as their gas and dust is exhausted (/ < 0.05). 
Today they would lie in the faint end of the DMF, mostly below the 
limits to which we can currently probe. They would need to be large 
stellar mass objects (since their stellar masses are already large at 
z — 0.5) but have little gas and dust today. They could plausibly 
be intermediate mass (logM, — 10.5 — 11.5) early type galaxies 
(ETG) in the local Universe, although they would still be relatively 
young since they were forming stars actively at z = 0.4 — 0.6. Such 
depleted objects could have had much more dust in the past with ra- 
tios of > 4 for the closed box scenario and the model with mantle 
growth proportional to SFR. In fact, the dust content of such galax- 
ies in the past could be even higher since the build up of a hot X-ray 
ISM in ETG rapidly destroys any remaining dust (e.g. Jones et al. 
1994). This is an attractive solution as severe outflows are then no 
longer required to reproduce the strong dust mass evolution seen in 
the DMF. Such a scenario predicts a population of early type galax- 
ies with moderate dust content and moderate ages (< 5 — 6 Gyr) 
as the last remnants of their ISM is depleted and the dust gradually 
destroyed. H-ATLAS has in fact discovered some promising can- 
didates for this transitional phase which are discussed in detail by 
Rowlands et al. in prep. 

Although a closed model does not reproduce the complexity of 
dust and metal growth within galaxies, we note that this elementary 
model including mantle growth predicts the highest dust masses for 
galaxies with the same initial gas mass and SFRs. Inflows and out- 
flows of material simply reduce the dust fraction in the ISM. A full 
treatment of the build up of metals in galaxies from stars of differ- 
ent initial masses further compounds this since relaxing the instan- 
taneous approximation would produce less dust at earlier times (at 
larger values of /). The difficulties we have in producing the ob- 
served dust evolution with this elementary treatment are thus only 
going to be exacerbated once a more complex treatment is adopted 
and therefore our conclusions about the requirements for higher 
yields and condensation efficiencies are conservative. To address 
the issues above, in particular, the importance of the star forma- 
tion history and the role of the IMF, a more complex model of dust 
and chemical evolution is required which allows mantle growth, 
destruction and even the shape of the IMF to depend on the star 
formation rate of galaxies. This is beyond the scope of this paper 



and the reader is referred to Gomez et al. (in prep) for a more com- 
plete investigation of the origin and evolution of dust in galaxies. 

6.3 Final caveat 

There is one important way in which the observed dust masses 
could be over-estimated; through the dust mass absorption coef- 
ficient K. This normalises the amount of emission from dust to the 
mass of material present and is dependent on the optical properties 
and shapes of the dust grains (for a more thorough review of the 
literature see Alton et al. 2004). The value of k used here is based 
on that measured in the diffuse ISM of the Milky Way (Boulanger 
et al. 1996; Sodroski et al. 1997; Planck Collaboration 2011a) and 
also on nearby galaxies by James et al. (2002). This value is some 
70 percent higher than that predicted by some models of dust, in- 
cluding the silicate-graphite-PAH model of Li & Draine (2001), 
but lower than those measured in environments where dust may be 
aggregated, icy mantles or 'fluffy' (Matthis & Whiffen 1989; Os- 
senkopf & Henning 1994; Rrugel & Siebenmorgen 1994). Latest 
results from Planck (Planck Collaboration 2011a) do see a varia- 
tion in the dust emissivity with temperature which is expected if 
there is grain growth in the ISM. It is thus not inconceivable that 
K could be globally higher in galaxies with larger fractions of their 
ISM in states which lend themselves to the growth of grains, or 
where larger fractions of grains have a SNe origin, or are under- 
going destruction by shocks. For example Ossenkopf & Henning 
(1994) show that in only 10'^ years of grain evolution in dense en- 
vironments (10^ — 10* cm"'^) the dust emissivity can increase by 
a factor ~ 5 due to the freeze out of molecular ice mantles and 
coagulation. The same authors also show that changing the ratio of 
carbon to silicate dust can change the emissivity by ^ 40 percent. 
Such a change in global dust composition could reflect the time de- 
pendence of evolution of various dust sources (e.g. SN-II dominat- 
ing in early time) or metallicity changes favouring O or C-rich AGB 
phases. The mechanism for changing the fraction of the ISM in the 
densest phases conducive to mantle growth could be triggered star 
formation and feedback (e.g. following an interaction). The frac- 
tion of gas in dense clumps has been found to increase markedly 
in parts of GMCs which are affected by feedback from recently 
formed OB stars (Moore et al. 2007). Draine et al. (2007) find that 
for local SINGS galaxies there is no need to consider ice-mantles 
in the modelling of the dust emission, but similar modelling has 
not been attempted for higher redshift and more sub-mm luminous 
sources such as the H-ATLAS sources. 

A measurement of k at Herschel wavelengths (but for local 
normal galaxies) has been attempted by Weibe et al. (2009) and 
Bales et al. (2010b). Both works, however, suggest a much lower 
value for k, which would increase the dust masses estimated here 
by a factor ~ 3. Given the already difficult task in modelling the 
dust masses, we do not believe that «;25o can be significantly lower 
than the values assumed here. A determination of k for H-ATLAS 
galaxies is ideally required (as these are sub-mm selected sources 
which may preferentially have higher k). Should an enhanced k at 
higher redshifts be the explanation for the large sub-mm luminosi- 
ties of H-ATLAS galaxies then this has important implications for 
the interpretation of all high-z SMG and Herschel observations. A 
change in k will lead to a change in the opacity of galaxies since the 
interaction of the grains with optical/UV photons will be altered. A 
strong test is to look at the effects of different k on the attenuation- 
inclination relation in the optical as differing values of n in the 
sub-mm will (for a fixed observed sub-mm flux) produce different 
values for the dust opacity in the optical-UV (see Popescu et al. 
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2011 for further details). For galaxies in the Millennium Survey 
(Driver et al. 2007) the Li & Draine (2001) values of k (which are 
lower than those used here by 70 percent) gave the best consistency 
with the observed attenuation-inclination relation, however it will 
be interesting to see the results of similar modelling for H- ATLAS 
sub-mm selected sources (Andrae et al. in prep). One result of an 
increasing k with redshift would be a flattening of the attenuation- 
inclination relation with redshift. 

A thorough investigation of all the implications using radia- 
tive transfer modelling is required but a change in k is likely to 
affect dust masses and the outputs of semi-analytic models which 
try to predict the SMG populations. If the FIR luminosity of high-z 
galaxies is not dominated by obscured star formation (i.e. there is 
a contribution from low opacity diffuse ISM or 'leaky' star form- 
ing regions) then a change in k may also lead to a bias in SFR 
estimated via FIR luminosities. Very high dust masses and sub-mm 
fluxes for SMG in the early Universe have proved challenging for 
dust formation models and semi-analytic models of galaxy forma- 
tion. In addition to exploring additional sources of dust and IMF 
variations to explain the SMG populations, it is worth considering 
of the possibility of dust grain property evolution as well. 



7 CONCLUSIONS 

We have estimated the dust mass function for the Science Demon- 
stration Phase data from the //er.sc7je/-ATLAS survey, and inves- 
tigated the evolution of the dust mass in galaxies over the past 5 
billion years. We find that: 

• There is no evidence for evolution of dust temperature out to 
z = 0.5 in this 250/im selected sample. 

• The dust mass function and dust mass density shows strong 
evolution out to z = 0.4 — 0.5. In terms of pure mass evolution this 
corresponds to a factor 4—5 increase in the dust masses of the most 
massive galaxies over the past 5 billion years 

• Similar strong evolution is found in the ratio of dust-to-stellar 
mass and V-band optical depth - //er.sc-/;e/-selected galaxies were 
more dusty and more obscured at z = 0.4 compared to today. 

• In order to account for the evolution of the dust content we 
need to radically alter chemical and dust evolution models. We can- 
not reproduce these trends with Milky Way metal or dust yields or 
star formation efficiencies. 

• H-ATLAS 250/im selected sources are highly efficient at con- 
verting metals into dust, either through mantle growth or through 
a bias in the IMF towards higher mass stars. They must also be 
observed following an episode of star formation (either recent for- 
mation or recent major burst) where the gas has been consumed at 
a much faster rate than galaxies like the Milky Way today. 

• As dust and gas (particularly molecular gas associated with 
SF) are tightly correlated in galaxies, this increase in dust content is 
suggestive of galaxies being more gas rich at 2: = 0.5. According to 
the simple chemical model, we are possibly witnessing the period 
of growth toward peak dust mass when gas fractions are ^ 0.5 
or higher. This strong decline in gas and dust content may be an 
explanation for the decrease in star-formation rate density in recent 
times as measured in many multi-wavelength surveys. 

This study uses only 3 percent of the area of the H-ATLAS 
data. Future improvements will come from the wider area coverage 
of the full survey, reducing uncertainties due to cosmic variance and 
small number statistics. Use of deeper optical/IR data from forth- 
coming surveys such as VISTA- VIKING, pan-STARRS, DES and 



VST-KIDS will also allow us to push to earlier times and higher 
redshifts to find the epoch of maximum dust content in the Uni- 
verse. 
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APPENDIX A: CHEMICAL EVOLUTION MODELLING 

This simple chemical evolution model describes the star, gas, metal 
and dust content of a galaxy making the instantaneous recycling ap- 
proximation. The mass fraction of metals, Z in this model changes 
as a mass ds of the ISM is formed into stars assuming no inflows 
or outflows via the following equation (Edmunds 2001): 

d(Zg) = apds + (1 - a)Zds - Zds (9) 

where g is the gas mass and a (Eq.llOt is the fraction of mass 
from a generation of star formation which is locked up in long-lived 
stars or remnants mjj as determined by the initial mass function 



[m — rnnirn)] <f){m)dm 



(10) 



p is the effective yield of heavy elements from stars p = 
p' /a 0.01 where a ~ 0.7 in agreement with Milky Way val- 
ues for a Scalo IMF. 

In a closed box model (i.e. no inflow or outflow of material), 
the total mass of the system (Mtot = gas + stars) is unity so 
that the fraction of gas in a galaxy (the ratio of gas mass to total 
baryonic mass) is / — g. In this scenario, the initial conditions are: 
Z = Q?ig = f = l and the gas mass of the galaxy is given by 
g = 1 — cts. The analytic solution for the metal mass fraction Z is 
(Eq.[n): 



Z = -pin/. 



(11) 



An early episode of star-formation prior to the evolution of the 
closed box would pre-enrich the gas and increase the interstellar 
metallicity (pre-enrichment is often invoked to explain the metal- 
licities of globular clusters in the Milky Way). We can include pre- 
enrichment of the ISM with metals Zi using 



Z = Z,- pin/ 



(12) 



where Zi 0.1 - 0.2 Zq (VanDalfsen & Harris 2004). Cor- 
respondingly, the dust mass fraction y varies with ds via: 



d{yg) = apxids + (1 - a)x2Zds - yds 



(13) 



where x is a parameter to describe the fraction of the mass of 
interstellar metals in dust grains from supernovae remnants or their 
massive star progenitors (xi)> and/or from the stellar atmospheres 
of low-intermediate mass stars (LIMS: X2)- The analytic solution 
is given in Eq. [T4]for y = at g = 1 and for a = 0.7 (typical 
locked up fraction for a Scalo IMF): 



y = 2.3 



(Xi-X2)(l -/"■«) 



+ 0.43x2 



pln(l/f) (14) 



ln(l//) 

For the special case where Xi = X2 = X' Eq.[T4]reduces to: 



y = xpin(i/f). 



(15) 



We can add an additional term to the dust mass from stars by 
assuming that grains accrete at a rate proportional to the available 
metals and dust cores in dense interstellar clouds (following Ed- 
munds 2001): 



y = XPln(l/ff) + e??c(2 - y) 



(16) 



where e is the fraction of the ISM dense enough for mantle 
growth (here we set this arbitrarily to 0.3), 77c is the efficiency of 
interstellar depletion in the dense cloud (i.e. if all the metals in the 
dense clouds are accreted onto dust grains then 77c = 1). 

Dust destruction via supernova shocks can be added to this el- 
ementary model by assuming some fraction S of interstellar grains 
are removed from the ISM as a mass ds is forming stars (therefore 
adding a term ~Sds to Eq.|13b. In this work, we test both a constant 
fraction with 5 = 0.3 (appropriate for MW-type galaxies and there- 
fore provides a testcase with a minimum destruction level expected 
for the H-ATLAS spirals) and a function that varies proportionally 
to the SFR (since a higher SFR equates to a higher Type II SN rate). 



Outflow 

We include a simple model for outflow of gas, in which gas is 
added or lost from the system at rates proportional to the star forma- 
tion rate. For large galaxies this outflow rate is assumed to be less 
than four times the SFR (A/a ^ 4; see Eales & Edmunds 1996 
for discussion; this corresponds to a galaxy which retains only ~ 
20 per cent of its original mass). We do not consider inflow of unen- 
riched material since this only slightly reduces the dust mass w.r.t. 
the closed box model and doesn't significantly change the evolution 
of a galaxy (Edmunds 2001). One can imagine a scenario with in- 
flow of pre-enriched material (e.g. merger), providing new material 
for star formation, even at later times when the original gas mass 
of the galaxy has been consumed through the star formation effi- 
ciency parameter k. Modelling the effects of this on the dust mass 
is beyond the scope of the simple model presented here. 

Outflows remove dust from the interstellar medium via 
— Xyds. The solution is given by Eq.[T7]if destruction 5 = 0: 



j/(outflow) = 



y 



(17) 



1 + X/a 

The gas mass g is related to the gas fraction / in this model 



by: 



<; (outflow) = 



/ 



l + (A/a)(l-/)' 
the metallicity mass fraction, Z: 



(18) 



Z(outflow) = -J^i^, (19) 



1 + A/a' 
and the total mass of the system is: 
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Mtot (outflow) = l+lhMl, (20) 
1 + A/a 

Dust and Stellar Mass 

The dust mass per unit stellar mass for the elementary model for 
equal x with no mantles, destruction or outflow, is given by Eg. 1211 

as 1 — <? 

We can rewrite Eq. [2T|as a function of time, since SFR ^{t) 
is related to the gas mass via is related to the gas mass via 

V'(t) = (22) 

where k is the star formation efficiency measured in Gyr~^ 
and the variation of g with time is 

High values of k will result in a higher SFR and a more rapid 
build up of the final stellar mass for the same initial gas mass. 

For outflow models, the dust mass fraction and the gas mass 
is reduced as described in Eqs. ?? - 1191 



